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ABSTRACT 

Overdensities in the distribution of low latitude, 2MASS giant stars are revealed by systematically peeling 
away from sky maps the bulk of the giant stars conforming to "isotropic" density laws generally accounting for 
known Milky Way components. This procedure, combined with a higher resolution treatment of the sky density 
of both giants and dust allows us to probe to lower Galactic latitudes than previous 2MASS giant star studies. 
While the results show the swath of excess giants previously associated with the Monoceros ring system in the 
second and third Galactic quadrants at distances of 6 - 20 kpc, we also find a several times larger overdensity 
of giants in the same distance range concentrated in the direction of the ancient constellation Argo. Isodensity 
contours of the large structure suggest that it is highly elongated and inclined by about 3° to the disk, although 
details of the structure — including the actual location of highest density, overall extent, true shape — and its 
origin, remain unknown because only a fraction of it lies outside highly dust-obscured, low latitude regions. 
Nevertheless, our results suggest that the 2MASS M giant overdensity previously claimed to represent the core 
of a dwarf galaxy in Canis Major (/ ~ 240°) is an artifact of a dust extinction window opening to the overall 
density rise to the more significant Argo structure centered at larger longitude (/ ^ 290° ±\Q°,b ^ -4° ± 2°). 

Subject headings: Galaxy: structure - Galaxy: disk - Local Group - galaxies: interactions 



L INTRODUCTION 

Newbere et al. '(•^OOB identified a low-/? excess of stars in 
the Sloan Digital Sky Survey (SDSS) at [IM - [223, +20]°, 
with distances beyond, and with a larger apparent scale- 
height than, the nominal Milky Way (MW) disk. They at- 
tributed this excess to the presence of either a newly dis- 
covered dwarf galaxy or tidal stream lying just outside the 
Galactic disk. Majewski et al. (2003, M03 hereafter) identi- 
fied the sa me structure in their study of 2MASS M giants. 
Ilbata et all ( ii2003 ) observed the structure with optical CMDs 
for fields spanning over 100° of longitude, and proposed as 
a most likely explanation of this apparent "ring" around the 
MW disk that it is "a perturbation of the disk, possibly the re- 
sult of ancient warps". Meanwhile, Yannv et al. (2003) relied 
on SDSS spectra and imaging in other areas to argue the struc- 
ture repr esents the rem ains of a tidally disrupted MW satellite, 
whereas fHelmi et all ( >200 3) explored models of the structure 
as either a "transitory localized radial density enhancement" 
from particles stripped off a satellite on a recent peri-Galactic 
passage, or the analogue of "shells" found around elliptical 
galaxies deriving from minor mergers several Gyr in the past. 

More detailed study of the distribution of 2MASS M gi- 
ants has shown that (1) this newly-found structure is in- 
deed ring-like, spanning > 170° around the disk at a typi- 
cal Galactocentric distance of ~ 16 kpc and broadened into 
both north and south Ga lactic hemisp heres ( Rocha-Pinto et al. 
l2003L lPaper I'hereafter:' Martin et al.k200 4a. M04a hereafter), 
(2) the structure contains stars enriched > lOx highe r tha n 
the [Fe/H]= -1.6 previously reported by Yannv et al. ( 2003b, 
belying a large metallicity spread ( Paper I: Crane et al. 2003h 
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like that seen in MW satellite galaxies, and (3) the stars in 
the structure follow a slightly non-circular orbi t with a rel- 
atively low (13- 20 km s~' ) velocity dispersion llCrane et alJ 
l2003t [Martin et alJ l2004bl) . Both globular and old open 
clusters spatially and kinernatically correlate d to the struc- 
ture have b een identifi ed (Crane et al. 2003; Frinchabov et alJ 
■2004; .Bellazzini et al.. 2004): these clusters follow a well- 
defined age-metallicity relation. This combined evidence sup- 
ports the view that the anticenter "ringt't' represents tidal de- 
bris from a Sagittarius (Sgr) dwarf-like galaxy along an orbit 
with low inclination to, and orbital radius just larger than, the 
MW disk. 

More recently, probing asymmetries i n the 2 MASS M gi- 
ant distribution about the Galactic plane, "MO^al claim to have 
found the core of a satellite galaxy in Canis Major ("CMa"; 
I ^ 240°) , which they arg ue to be the progenitor of the ring. 
However, iMomanv et alJ (12004) have argued that CMa may 
be the MW disk stellar warp, based on comparison of ob- 
served to synthetic CMDs generated by models of a warped 
MW disk. This has prompted vigorous d ebate over the 
nature of C Ma (Bellazzini et al. 2004; Martin et al. 20041^ 
ISbordoneet al. 2005; Penarrubia et al. 2005). Whatever its 
origin, that there is a "CMa overdensity" has apparently 
gained wide acceptance (e.g., Forbes et al. 2004' 'Tosi et ajj 
2004; Sbordone et al. 2005; Martinez-Delgado et al. 2005). 

In this contribution we call into question a "core" overden- 
sity in CMa, and find an apparently more significant giant star 
overdensity at / ^ 290°. 

2. DATA AND COMPARISON MODEL 

We extract a giant star candidate sample from the 2MASS 
all-sky point source release as all stars having dereddened 
0.85 < J-Ks < 1 .5, Ks < 13.0, photometric quality flag AAA, 
and that meet the M giant color locus restriction used by 
M03. For each star we derive a distance probability distribu- 
tion function (DPDF) based on its dereddened, PSF-fitted Ks 
and J-Ks, the va riation of the color- absolut e magnitude rela- 
tion with [Fe/H] dlvanov & Borissoval2002^ . and an assumed 
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[Fe/H] distribution for giants typical of the ring (a Gaussian 
with mean [Fe/H]= -1 and spread 0.4 dex). C alculation of the 
DPDF is described in more detail in iPaoer 1 We adopt as a 
single distance estimate, r, for each star the mode of its DPDF. 
Extinctions were calculated as {Aj,AkstEj-Ks) = 

f.9Q,Q.25,Q.65)aEB-v, where Eb-v is from 
:hle2el, Finkbeiner, & Davis (1998) and a = | 
Locha-Pinto et al. 2004). As with previous, similar 
studies (Paper 1; M04a) our analysis is limited to Eb-v < 0.55 
to avoid false overdensity signals from small-scale differ- 
ential reddening or errors in its estimation. However, we 
make several improvements over previous work: ( 1 ) To more 
closely follow the Eb-v = 0.55 "border" our maps are binned 
into finer, 1° x 1° (l, b) bins t han the 4° x 2° bins of M04a 
or the 4° x 4° bins of 'Paper f. (2) 'M04a depended on direct 
comparisons of one region of the sky to another presumed to 
be "symmetric" — a procedure that requires both regions to 
have Eb-v < 0.55 and that ultimately removes large regions 
of "good" sky from consideration because comparison fields 
are too reddened. Here we rely on a Galactic model of 
the number of stars from canonical Galactic populations 
as a function of distance and il,b) as a reference for all 
Eb-v < 0.55 fields, which allows us to probe areas of the sky 
previously ignored. 

The densities of the asymmetries sought are small com- 
pared to the MW foreground density. The goal of creating a 
Galactic model is to "fairly" remove the majority of this fore- 
ground to highlight the residual. Rocha-Pinto et al. (2006; 
in preparation) describes the Galaxy model in detail; here we 
outline the general procedures used to construct it. The key to 
such a model is that it be smooth, have a cylindrically sym- 
metric density distribution about the Galactic center, and be 
constrained globally by the observed density distribution in 
available, low Eb-v lines of sight, so that, in the mean, most 
of the 2MASS giants are removed when the model is sub- 
tracted from the "observed" {l,b,r) distribution to leave be- 
hind a high contrast picture of "overdensities". To this end, 
the model does not need to be exact, nor should it necessar- 
ily be interpreted as a true description of the MW giant star 
distribution, though we start from a parameterization meant 
to produce a classical "starcount model" from thin disk, thick 
disk and halo density laws (those of Sieael et al. 2002). These 
density laws are multiplied by adopted luminosity functions 
for the corresponding stellar populations and integrated from 
the brightest to faintest absolute magnitudes of stars that can 
be in our observed sample. After integration over solid an- 
gle (1° x 1°) and normalization of the results to the observed 
stellar density at longitudes not contaminated by known stel- 
lar streams or the Magellanic Clouds, we have the expected 
number of stars coming from each contributing MW popula- 
tion within a given r range. 

Because the 2MASS stars have distances calculated as if 
they have come from a hypothetical stellar population with 
([Fe/H]) = -1.0 and (7[Fe/H]= 0.4, stars from the simulated 
MW have their distances recalculated with a DPDF in the 
same way — i.e. as if their simulated metallicity were un- 
known. That our model MW succeeds in removing the bulk 
of the canonical MW populations of giants is shown in the dif- 
ference between the observed and model 3D density distribu- 
tions (Fig.[0, where stellar densities in the projected distance 
ranges r < 6 kpc and 20 < r < 50 kpc are essentially elim- 
inated. In contrast, in the distance range dominated by the 
ring, 6 < r < 20 kpc, significant residual densities — which 
we claim are true overdensities — can be seen. 
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Fig. 1 . — Residual 2MASS late type giant star density with respect to our 
MW model projected onto the celestial sphere at different distance ranges. 
Contour levels represent geometrically increasing densities in number of stars 
deg"-. The solid lines show Eb-v contour levels at 0.55 and 1.20 mag; re- 
gions having Eb-v > 0.55 have been discarded from the analysis. This figure 
supports our claim that the "CMaf t' overdensity is only an outlying part of 
the larger Argo overdensity seen in the 6 < r < 20 kpc range. 



3. THE ARGO OVERDENSITY 

Most of the "^"-shaped arc of overdensity sweeping all 
of the way across the middle panel of Figure 1 can be at- 
tributed to the "anticenter ring/warp/Monoceros stellar struc- 
ture," discussed by the various studies referenced in § 1 . Here 
we focus specifically on the very obvious, large overdensity 
in the third and fourth quadrants from about Z = 210° to at 
least the left edge (/ = 310°) of our sample range and peak- 
ing at [l,b] ^ [290,-4]°. From the overall rise in density 
towards this specific point in Figure 1, we surmise that the 
trend continues and an even higher overdensity peak may lie 
somewhere in the obscured region with -4° < b < +2° near 
285° < I < 300°, which spans the constellations Carina, Vela 
and Puppis. Because the peak of this large overdensity re- 
mains uncertain and possibly obscured, we call the structure 
"Argo," after the large, ancient constellation later dismem- 
bered into Carina, Vela and Puppis by the lAU. As may be 
seen, the density of this region is significantly higher than 
any other visible giant star excesses — e.g., it is 4 to 8 times 
more dense than the densest, visible parts of the Northern and 
Southern Ring in the second quadrant. 

The overall impression is that this large excess of giant stars 
constitutes one coherent structure with a "core" at I ~ 290° 
but spreading well beyond. This core appears to be local- 
ized in distance. Figure |2] shows a 2MASS Hess diagram of 
the Argo core, created as the difference in Hess diagrams for 
fields at [l,b] = [285,-5]° and [l,b] = [285,-1-5]°. From this 
difference a well-defined red giant branch (RGB) associated 
with Argo emerges from the "RGB-smear" created by closer, 
foreground MW disk stars. The RGB residue has a highest 
Hess di agram density that is fit by an [Fe/H]= -0.7 population 
dlvanov & Borissova 200 2) at r ^ 13.8 kpc, but isochrones 
with a spread of metallicity not unlike that previously sug- 
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Fig. 2. — 2MASS CMDs for (top) a field with an excess of stars associated 
witli the Argo overdensity, (middle) a control field (same /, opposite b), and 
(bottom) their difference. The line in the top panel shows that the bulk of 
MW giants can be fitted by an [Fe/H]= -0. 1 RGB isochrone having a median 
distance modulus of 14 mag (i.e. 6.3 kpc). This isochrone is reproduced in 
the lower panel as a dashed line to demonstrate that the residual RGB we 
associate with Argo is mostly fainter. The soHd lines in the bottom panel 
show [Fe/H]= -1.4, -0.7 and +0.0 RGB isochrones placed 13.8 kpc from the 
Sun to match the distinct RGB from Argo as well as the previously suggested 
(see §1) large [Fe/H] spread of the anticenter ring. Residue al J-Ks ~ 0.5 
may represent Argo red clump stars. 



gested for stars in the ring (see §1) and at this same distance 
seem to reasonably account for the triangular spread of resid- 
ual red stars. That a red clump predominantly at fainter mag- 
nitudes remains after Hess diagram differencing is not incon- 
sistent with the presence of an Argo red clump at the same 
general distance dSalaris & Gir ardi 2002). 

Another aspect of interest is that the isodensity contours 
tend to have ellipsoidal - not warp-like — character, with the 
major axes of the isopleths "jutting" out obliquely from the 
MW disk. In this respect, Argo resembles a large, distorted 
dSph galaxy, similar to the Sgr system. If, like Sgr, Argo is a 
coherent, symmetrical structure, we can "fitt't' ellipses to the 
visible isopleths to hypothesize its overall appearance (Fig. 
3). The "eyeballt't' fit ellipses are centered at one point in the 
Eg-v > 0.55 region, but have variable axis ratio and inclina- 
tion angle. The suggestion of a "twisting isophote" character 
resembles the M giant distribution of the Sgr core (see, e.g., 
M03's Figs. 7e and 7f), which, with its own "ring" of tidal de- 
bris wrapping around the MW, may serve as a useful paradigm 
for Argo. The apparent inclination of Argo to the Monoceros 
ring (which seems to have a closer alignment to the MW disk) 
recalls the fact that the major axis of the Sgr core is canted (by 
6°) with respect to its tidal stream (M03). These structural 
properties of the Argo system are certainly tantalizing support 
for the notion that it may be a tidally disrupting dSph-system 
— perhaps the progenitor of the Monoceros ring. 

In theu- 2MASS M giant studv. -M04a also identified an ap- 
parently ellipsoidal overdensity, but with a geometrical center 
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Fig. 3. — Same as middle panel of Fig. 1, but with greyscales fit by the el- 
lipses shown. The dashed line indicates the Eg-v = 0.55 reddening boundary 
adopted in the M04a study. The dotted and solid lines show the Eg-v = 0.55 
and 1.2 bound aries, respectively, from a finer extinction map resolution than 
used bv iM04aL 



at {l,b] ^ [240,-8]°. From this core-like feature, at an appar- 
ent distance of 10.0 < Rcc < 15.8, M04a concluded that the 
center of the Monoceros system is in CMa. In contrast, our 
map shows nothing particularly special about this region, ex- 
cept that it lies within the overall density rise towards the Argo 
core, which in our map is 2-4x more dense t han the excess 
at the position of the CMa "core." Since the lM04al analy- 
sis is using ostensibly the same database as ours, it is critical 
to understand the source of the discrepancy in results. Clues 
may lie in differences in how reddening affects each analy- 
sis: iM04a| used both a coarser Eg-v map resolution and were 
forced to discard regions with highly extinguished Galactic 
plane-mirrored counterparts that we preserve thr ough u se of 
our model. The decreased area available to the M04a' anal- 
ysis is demonstrated by the difference between their and our 
Eb-v = 0.55 boundaries (Fig. 3). Much of Argo, including the 
highest density patch , has be en removed from their analysis. 
In addition, with the lM04al boundary it becomes less obvi- 
ous that the overdensity at I ^ 240° is related to the one at 
I ~ 290° because of the intrusion of the elongated dust feature 
at I ~ 265° extending south to ~ -12°. The 'aproximately 
ellipsoidal' shape of CMa observed by these authors may be 
an artifact caused by the neighboring very reddened map cells 
— the 265° reddening "fingert'f and another at I ~ 220°; the 
likely coordinates for the Monoceros structure core as given 
by M04a look to be the geometrical center of the outer Argo 
overdensity seen through this reddening window in their map. 
Since the same arguments might be made about the nature and 
likely center of the Argo overdensity derived from our maps, 
it is important to stress both the provisional nature of our as- 
sessment of Argo and that most important clues about this star 
system still remain hidden. 

4. CONCLUDING REMARKS 

Our analysis of the 2MASS giant star distribution in the 
outer Galaxy, gives strong evidence for the presence of an 
excess of stars, in the distance range of the Monoceros ring, 
stretching across at least 210° < / < 310° with a concentration 
in the Argo region {I ~ 290°), and suggests the CMa "core" is 
a reddening artifact sitting on the stern of Argo. At minimum, 
our results from an improved analysis of 2MASS M giants in 
high Eb-v regions show a different structure to the stellar den- 
sity enhancements in the MW third and fourth quadrants than 
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Fig. 4. — Comparison between the gas warp and the stellar overdensities 
found in this paper. The dai'k and light contour Hnes (blue and red lines, in 
the electronic version) in this plot refers to southern and northern hemisphere 
asymmetries, respectively, projected onto the Galactic plane. Contour levels 
indicate 0.2, 0.4, 0.6 and 0.85 of the maximum gas (5.36Mq pc"^) or stellar 
(4554 M giants kpc"' rad"') cell density in this plot. The position of the Sgr 
dSph, as well as CMa and Argo, is shown. 



previously reported. Figure 3 further explores our discovered 
peak overdensity in Arg o as if it (not CMa) were the core of 
an elongated, possibly disrupting dwarf satellite galaxy, per- 
haps the source of the Monoceros ring. But further work is 
necessary to verify this and alternative explanations. For ex- 
ample , based on the distribut i on of 2MASS K giant candi- 
dates, iLopez-Corredoira et alJ l|2003) report the discovery of 
a stellar warp in the outer disk. They base this conclusion on 
an apparently sinusoidal excess of stars with respect to /, un- 
der equatorial symmetry. The south-north asymmetry in their 
data is strongest at I ^ 245°, coinciding with CMa, but the 
model proposed by the authors to reproduce these data peak 
around / ^ 290°, where Argo is found. While this may give 
the impression that Argo is the stellar component of the MW 
warp, their analysis relies on the same symmetry comparison 
methodology that affected other previous (e.g.. Paper I, M04a) 
low latitude 2MASS studies and, as such, cannot properly ac- 
count for the presence of any satellite near and/or spanning 
across the MW mid-plane. 



It is also useful to compare the HI Galactic warp with the 
stellar asymmetries to see whether these coincide. Investiga- 
tions of the gas warp show strong deviations from the Galac- 
tic plane both near / ~ 245° and / - 290° (see Fig. 5 by 
ISodroski et al. 1987). Fi gure|4]compares the stellar asymme- 
tries found by this work with the asymmetries in the gas warp, 
according to the 3D data by Nakanishi & Sofue (2003). Al- 
though the southern gas warp asymmetry is similar in shape 
to that coming from the stars in the third and fourth Galac- 
tic quadrants, the gas warp is 2 and 5 kpc farther from the 
Sun than the stellar overdensity at the longitudes of CMa and 
Argo, respectively. Moreover, quite contrary to what is seen 
in the stellar overdensity, the HI warp shows a symmetrical 
deviation to the Northern hemisphere across the Z = - 180° 
Galactic meridian and, moreover, it is densest in this north- 
ern hemisphere. Since there is no counterpart to the north- 
ern hemisphere gas warp in the stellar data, we conclude that 
Argo is not likely to be simply a stellar counterpart to the 
gas warp. On the other hand, it is possible that the pres- 
ence of this putative satellite galaxy could be the cause of the 
HI warp, given that interactions with satellite dwarf ga laxies 
is the leading explanation for disk warp s ( Shang et al. 119981: 
ISchwarzkopf & DettmaJ200lHReshetnikov et al..2002h . 
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